Asymptotic behavior of the Kleinberg model 
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We study Kleinberg navigation (the search of a target in a d-dimensional lattice, where each site 
is connected to one other random site at distance r, with probability ~ r~ a ) by means of an exact 
master equation for the process. We show that the asymptotic scaling behavior for the delivery 
time T to a target at distance L scales as T ~ In 2 L when a — d, and otherwise as T ~ L x , with 
x — (d — a)/(d + l — a) for a < d, x — a — d for d < a < d + 1, and x = 1 for a > d + 1 . These values 
of x exceed the rigorous lower-bounds established by Kleinberg. We also address the situation where 
there is a finite probability for the message to get lost along its way and find short delivery times 
(conditioned upon arrival) for a wide range of a's. 
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In a now famous study [l( the social psychologist Stan- 
ley Milgram asked randomly chosen people to send a 
postcard to a disclosed target in the USA. The partic- 
ipants were to send the card to a person they knew on 
a first-name basis, who will then send it on to another 
acquaintance, and so on, until it reached the desired des- 
tination. 20% of the cards successfully reached the target 
using, on average, chains of 6.5 acquaintances, confirm- 
ing the notion that the network of social contacts has the 
small-world property Q: a very short path, typically log- 
arithmic in the size of the system, connects between any 
two nodes. How does the message find its way, let alone 
so efficiently? Searching and navigation problems such 
as this 0, |, i, i, 0, S II 03, 03, El G3 are relevant to 
several disciplines, from sociology, to efficient algorithms 
in computer science, to the understanding of the foraging 
behavior of insects and animals. 

The problem has been elegantly addressed in the sem- 
inal work of Jon Kleinberg 0, 0]. Kleinberg considers 
an (L x L)-square lattice, where in addition to the links 
between nearest neighbors each node i is connected to 
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a random node j with a probability pij 
(rij — \ri— Tj\ is the Euclidean distance between nodes 
and j, and the sum in the denominator excludes k = i). 
Suppose that a message is to be passed from a "source" 
node s to a "target" node t, along the links of the net- 
work, by a decentralized algorithm — an algorithm that 
relies solely on local information. Kleinberg shows that 
when the exponent a equals d, the lattice dimensionality, 
an algorithm exists that requires less than In 2 L steps to 
complete the task. If a d, the delivery time, T, grows 
as L x , with rigorous lower bounds, x > xk 0,0,03; 
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< a < d, 
a > d . 
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Moreover, no local algorithm performs better, function- 
ally, than the simple-minded greedy algorithm: Pass the 
message forward to the neighbor node that is closest to 
the target (geographically). 

In this letter we study the asymptotic long-time be- 
havior of the Kleinberg search process. We find the exact 
form of the exponent x(a), and we show that T ~ In 2 L is 
the actual scaling (not just a bound) for the special case 
of a = d. Our approach is based on a master equation for 
the full probability distribution for completing a search 
within a given time. This formalism also enables one to 
consider the possibility of the message getting lost along 
its way, and we discuss briefly some surprising outcomes 
of that scenario. 

Because the message gets closer to the target with 
each successive step, nodes are never revisited and one 
can view each long-range link as being created only as 
the message arrives at its site of origin. The message 
closes in on the target in a peculiar kind of directed Levy 
walk, consisting of a mix of "short" (one lattice spacing) 
and "long" (power-law distributed) steps. In this search 
the dimensionality of the lattice enters into consideration 
mainly as a Jacobian in the various sums (or integrals) of 
the equations describing the process. We therefore limit 
the following derivations to one dimension and general- 
ize the results for higher dimensions, having made the 
necessary adjustments. 

For convenience, to render the independent of i, we 
adopt periodic boundary conditions. Specifically, con- 
sider a ring of length 4L, with the source at and the 
target at L, and the range of the long-contact links lim- 
ited to 2L — 1 [lj| . In that case, 
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where A is a normalizing factor. Let P{n; I) denote the 
probability that a message, at distance I from the tar- 
get, takes n additional steps to reach the target. Once 
the message is at the target it takes no additional time 
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FIG. 1: (Color online) Mean delivery time, Tl, as a function 
of the long-contact exponent, a. Note the perfect agreement 
between simulations (solid line) and the results from Eq. Q 
(symbols). Shown are results for three values of L. Inset: 
Distribution of the delivery time for the case of a = 1 and 
L — 1000 as computed from (|3} (solid line) is compared to 
simulations (symbols). 



to reach it, so P(n; 0) = <5 ni o- Likewise, P(0;l) = <5o,f, 
since the only way to reach the target in steps is if the 
message is already there to begin with. 
P{n\ I) satisfies the equation 
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The first term on the rhs represents the events that the 
first of the additional n steps is a long step of length k, 
in which case the message would come to within distance 
\l—k\ from the target. The second term represents a short 
step, that advances the message a single lattice spacing. 

For our main purpose here it is sufficient to consider 
just the first moment (n) = Ti, that is, the mean deliv- 
ery time from a site a distance I away from the target. 
Multiplying Eq. ([3]) by n and summing over n, we get 
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for I = 1, 2, . . . , L. Numerical integration of Eqs. © and 
@ yields perfect agreement with the results from direct 
simulation of the Kleinberg navigation process on a ring 
(Fig. DO). 

Using the fact that T-k = and defining Dk — — 
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for i = 1, 2, . . . , L — 1, while for I = we have D\ — 
Ti - T = 1. 

As a quick check, consider the limit of a — > oo, when 
all the long-range contacts are restricted to length 1 , and 
therefore one expects T\=l. Indeed, in this case all the 
k~ a terms in the equation tend to zero, unless k = 1, and 
we get Di + i — Di — 0, which along with D\ = 1 yields 

Dk = 1, and Ti = J2k=i ^ k = ^ J us ^ as expected. 

Next, consider the opposite limit, of a — 0, where the 
distribution of long-rage contacts is homogeneous. In this 
case A = [2(2L — 1)] _1 and we obtain from ([5]), 
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Although this equation can be solved exactly, a continu- 
ous approximation, 
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assuming L > I ^> 1, works just as well. In view of 
the boundary condition D(0) — 1, this has the solution 

D(l) = exp(-/ 2 /4L). Then, T(L) = J Q L D{1) dl. The up- 
per integration limit may be safely replaced with oo, due 
to the rapid decay of the gaussian, and a simple change of 
variables yields T(L) ~ L 1 / 2 , in perfect agreement with 
the Kleinberg bound for a = 0. 

For larger values of a we are not as fortunate as to 
find a full analytic solution, but we can still obtain 
the asymptotic behavior. For < a < 1 we take a 
hint from the solution for a = and make the ansatz 
D(l) = f(l /L), where f(x) is a smoothly decreasing 
function; f(x) = 0(1) for x < 1, and decays very rapidly 
(e.g., exponentially) for x > 1. Consistent with this be- 
havior, the derivative at the crossover point x* = 1 is 
/'(a;*) = —0(1). This ansatz is nicely confirmed by nu- 
merical integration of Eq. (|5|) . 

Apply now Eq. $5$ to the crossover length Z» = L 1 ^. 
The lhs is 
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while the sums on the rhs can be estimated by replacing 
Di with a constant for I < Z„, and zero for I > Z», yielding 
-Al\- a . But A - 1/L 1 - 01 , leading to -1/(3 = (1 - 
a)//3 - (1 - a), or /3 = (2 - a)/(l - a). Finally, 
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' f(x)x 1 ^' 1 dx ~ L 1 ^ . 
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T(L) ~ i( 1_Q )/( 2_Q ) ( 0<a<l. 
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For a > 1, we sum Eq. |5]) over Z, taking into account 
that D\ = 1, and rearrange: 

l k 

D 1+1 -1 = -AJ2 (l + m)- a D k . (7) 

fe=l m=~k+2 

This can be obtained more directly also by rearrang- 
ing Eq ((3]). The inner sum over m can be approximated 
by an integral, yielding 
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Since a > 1, A converges, as L — > oo, and we may simply 
follow powers of I. Assume first that Di ~ I' 13 . The dom- 
inant — 1 on the lhs of the equation must be balanced by 
the dominant term on the rhs, which scales as —l 2 ~ a ~ 13 , 
so (3 — 2 — a. This leads to 
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The upper limit on a follows from the fact that Tl cannot 
increase faster than linearly in L. 

For the special case of a = 1, we have 
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and A ~ (In L) 1 . To counter the — 1 on the lhs one must 
then allow that D\ ~ (lni)/Z, for large I, which leads to 
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exactly as the Kleinberg upper limit for this case. 

The dominant —1 term on the lhs of ([8} might also be 
cancelled if D(l) = l—g(l), where g(l) vanishes as I — > oo. 
Substituting this ansatz in ([8]) we find g(l) ~ l 2 ~ a , and 
integration of D(l) yields 
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In this case, the condition a > 2 prevents the correction 
term from growing faster than linearly. 

In Fig. (2k- we compare our various results for x(a) in 
one dimension to computer simulations (l6( | and to the 
Kleinberg bounds, Xk (a)- For a — > 2 we expect logarith- 
mic terms, as the main contribution and the correction 
term in (jTTJ) approach then the same power-law. The 
logarithmic behavior makes it very difficult to extract re- 
liable estimates from simulations for the exponent x(a) 
near a = 1 and 2. 

The foregoing results, argued for one-dimensional lat- 
tices, are easily generalized to higher dimensions: the 
space dimensionality, d, enters in the various integrals 
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FIG. 2: (Color online) Delivery time exponent, x, as a func- 
tion of a in (a) d = 1 and (b) d — 2 dimensions. Results from 
Eq. (|12|l (solid lines) are compared to the Kleinberg bounds 
(broken lines) and to simulations (symbols). 



(or sums) through the Jacobian for rf-dimensional inte- 
gration. Making the necessary adjustments, we find 
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Once again, these results agree with the Kleinberg 
bounds x > xk- Eq. (|T2|) fits simulations results for 
d = 2,3,4 nicely (we did not test higher dimensions). 
Results for d = 2 are shown in Fig. [3t>- 

Returning now to the probability distribution, P(n; I), 
a standard way to tackle Eq. © is through the generating 
function 
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The "fugacity" z may be interpreted as the probability to 
complete a single step successfully. In that case P(z;l) is 
the total probability to complete the delivery successfully 
(to a target at distance I). We defer a more detailed study 
of P(n] I) for future work and focus for the moment on 
the conditional average of the delivery time, T z ^i, that 
is, the average delivery time conditioned upon successful 
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FIG. 3: (Color online) Conditional delivery time, T z> i, as a 
function of a, for L — 200 (squares), 400 (circles), and 600 
(triangles). Results for z = 0.9 (empty symbols) are compared 
to perfect transmission (solid symbols). 



arrival at the target: 

Y, n = Q riP{n;l)z n _ dP(z;l)/dz 
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Multypling ^ by z n and summing over n, minding 
the boundary conditions, we obtain 
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Taking the derivative of this equation with respect to z, 
and writing dP(z; l)/dz = F(z; I), we have 
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F(z; I) = A ^ k- a [zF(z; \l - k\) + P(z; \l - fc|)] 
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(l - A ^ k ~ a ) l zF ( z '> 1 ~ !) + p ( z '> 1 - l ) 



We are now ready to address the conditional average, 
at least numerically. Fixing the value of z, we solve 
Eqs. (fl4|) and (fT5|) recursively, up to ( = L, and use 
Eq. (fT5)l to compute the average. Typical results are 
shown in Fig. [3l For a fixed distance L the delivery time 
is remarkably small throughout a wide range of long- 
contact exponents, a < oj»(z;L), and saturates rapidly, 
T z .l ~ L, as soon as a exceeds a*. Thus, it seems that 
in the presence of losses small-world behavior does not 
imply a particular value of the long-contact exponent. 
If anything, there appears to be a small local maximum 
(barely perceptible in the figure) around a value of a that 
creeps towards d as L increases. These results seem rele- 
vant to the Milgram experiment, where it was estimated 
that z « 0.75 [ill. 

In summary, we have found the actual delivery time 
exponent in the Klcinberg navigation model, for the non 
small-world cases of a ^ d, and we have confirmed that 
T ~ In 2 L is the actual scaling (not just a bound) for the 
special case of a = d. We have also introduced master 
equations that allow one to study the Kleinberg model 
analytically in greater generality, and have looked briefly 
into the delivery time in the case of imperfect transmis- 
sion, when there is a finite probability for the message to 
get lost along its way. The distribution of delivery times 
and the analytical treatment of imperfect transmission 
are important problems left open to further inquiry. 
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